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Abstract 



We have investigated the spectral characteristics of a sample of bright 7-ray bursts de- 
' tected with the 7-ray burst sensors aboard the satellite Ginga. This instrument employed 

a proportional and scintillation counter to provide sensitivity to photons in the 2 - 400 keV 

0^ 



region and as such provided a unique opportunity to characterize the largely unexplored 
X-ray properties of 7-ray bursts. The photon spectra of the Ginga bursts are well described 



by a low energy slope, a bend energy, and a high energy slope. In the energy range where 
O ' they can be compared, this result is consistent with burst spectral analyses obtained from 

' the BATSE experiment aboard the Compton Observatory. However, below 20 keV we find 

' evidence for a positive spectral number index in approximately 40% of our burst sample, 

with some evidence for a strong rolloff at lower energies in a few events. There is a corre- 
^ ' lation (Pearson's r = -0.62) between the low energy slope and the bend energy. We find 

■ that the distribution of spectral bend energies extends below 10 keV. There has been some 

concern in cosmological models of GRBs that the bend energy covers only a small dynamic 
range. Our result extends the observed dynamic range and, since we observe bend energies 
down to the limit of our instrument, perhaps observations have not yet limited the range. 
The Ginga trigger range was virtually the same as BATSE's, yet we find a different range 
of fit parameters. One possible explanation might be that GRBs have two break energies, 
one often in the 50 to 500 keV range and the other near 5 keV. Both BATSE and Ginga 
fit with only a single break energy so BATSE tends to find breaks near the center of its 
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energy range and we tend to find breaks in our energy range. The observed ratio of energy 
emitted in the x-rays relative to the gamma-rays can be much larger than a few percent 
and, in fact, is sometimes larger than unity. The average for our 22 bursts is 24%. We also 
investigated spectral evolution in two bursts. In these events we find strong evidence for 
spectral softening as well as a correlation between photon intensity and spectral hardness. 
We also find that the X-ray signal below 30 keV itself softens in both of these events. There 
is one example of a strong x-ray excess at low energy. In addition to providing further 
constraints on 7-ray burst models, the description provided here of burst spectra down to 
2 keV should prove useful to future planned efforts to detect bursts at X-ray energies. 
Subject Headings: Gamma Rays: Bursts - X-Rays: General 
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1. Introduction 

Twenty-five years after their discovery, 7-ray bursts continue to defy explanation. Anal- 
ysis of burst energy spectra remains one of the principal methods for determining the phys- 
ical processes responsible for these events. Recent analysis of spectroscopy data from the 
Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observa- 
tory (CGRO) has demonstrated the diversity of burst spectral continua in the « 30 — 3000 
keV range (c/. Band et al. 1993; hereafter B93). However, few instruments to date have 
probed the X-ray regime of 7-ray bursts between 2 and 20 kcV. 

The first detections of X-rays in the 1-8 keV range from 7-ray bursts were reported by 
Wheaton et al. (1973), Trombka et al. (1974), and Metzger et al. (1974). Later results from 
the P78-1 (Laros et al. 1984) and Hakucho (Katoh et al. 1984) satellites confirmed that 
7-ray bursts often produce significant X-ray emission. In addition, these satellites provided 
the first evidence suggesting that the X-ray emission might outlast the main 7-ray event 
in some bursts (so called X-ray tails). These early data provided only modest spectral 
resolution in the X-ray regime, however, based on their results, the 7-ray burst detector 
(GBD) flown aboard the Ginga satellite was specifically designed to investigate burst spectra 
in the X-ray regime (Murakami et al. 1989). Ginga was launched in February of 1987, and 
the GBD was operational from March, 1987 until the reentry of the spacecraft in October, 
1991. Several important results have emerged from the study of burst spectra recorded with 
the Ginga GBD. Harmonically spaced line features have been reported from several Ginga 
events (Murakami et al. 1988; Fenimore et al. 1988; Yoshida et al. 1992; Graziani et al. 
1992), and have been interpreted as due to cyclotron absorption in the strong magnetic 
field of a neutron star. Recently, preliminary results indicate that BATSE has also detected 
absorption lines (Briggs et al, 1997). In addition to line features, Ginga has observed X-ray 
tails in a number of bursts, as well as x-ray preactivity in one event (Yoshida et al. 1989; 
Murakami et al. 1991; & Murakami et al. 1992b). More recently, the BcppoSax satellite has 
observed several bursts in both x-rays and gamma-rays (Piro, et al. 1997) and discovered 
soft x-ray afterglows (Costa et al., 1997). This latter discovery has opened the way for 
the long sought GRB counterparts, including one with a measured redshift (Metzger et al., 
1997). Here, we provide for 22 GRBs a compilation and description of 7-ray burst spectral 
continua emphasizing the X-ray spectral characteristics. 

Before BATSE, the analysis of 7-ray burst spectra suggested that the continua can be 
fitted adequately by a range of models, with power law, optically thin thermal bremsstrahlung 
and thermal synchrotron formulae providing acceptable fits to many spectra (c/. Mazets 
et al. 1982; Matz et al. 1985; and Hurley 1989). Although these functional forms ade- 
quately describe many burst spectra it is not at all clear that the corresponding physical 
processes are responsible for the burst spectra, because spectral shapes inconsistent with 
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each other may often be consistent with the same data. The compilation by B93 of av- 
erage spectra from a large sample of BATSE bursts provides the most detailed statistical 
description presently available of burst continua in the range from f« 30 — 3000 keV. (Also 
sec Schacfer et al. 1994.) It is our purpose here to provide an extension of this statisti- 
cal description of burst spectra down to 2 keV using an analysis of bright bursts observed 
by the Ginga GBD. In addition to providing further constraints on proposed burst source 
models it is our hope that the information provided here will be useful to future soft X-ray 
experiments which will require the range of GRB spectral characteristics to estimate, for 
example, the expected rate of burst detections. To facilitate comparison of our results with 
the pubhshed BATSE spectra we adopt the same burst spectral form as described in B93. 
Although it docs not imply any particular physical process, this form provides an adequate 
description of BATSE bursts and the Ginga bursts studied here. Given the excellent statis- 
tics generated by the BATSE detectors, the fact that the B93 shape can adequately fit the 
GRB spectrum is a strong argument that it should be used even though it does not imply 
any physical process. 

We begin in §2 with a brief description of the GBD instruments, as well as a discussion 
of the analysis techniques employed, including a discussion of the constraints placed on 
us by the lack of positional information for many of the Ginga bursts. In §3 we describe 

the spectral properties of the Ginga sample and compare these results with the previously 
reported BATSE spectra. In §4 we investigate the spectral evolution of a smaller number of 
bright bursts. In §5 we summarize the x-ray emission relative to the gamma-ray emission. 
We conclude in §6 with a summary and discussion of our principal results. 

2. Instrumental Summary and Data Analysis Techniques 

The GBD aboard Ginga consisted of a proportional counter (PC) sensitive to photons in 
the 2-25 keV range and a scintillation counter (SO) recording photons with energies between 
15-400 keV. Each detector had an !=a 60 cm^ effective area. The PC and SC provided 16 
and 32 channel spectra, respectively, over their indicated energy ranges. The detectors were 
uncollimated except for the presence of shielding to reduce backside illumination and the 
mechanical support for the PC window. The field of view was effectively vr steradian for 
both the PC and the SC. In burst mode the GBD recorded spectral data from the PC 
and SC at 0.5 s intervals for 16 seconds prior to the burst trigger time, and for 48 seconds 
after the trigger. This Memory Read Out (MRO) data was used for many of the spectral 
fits described here. In the event that MRO data was not available for a burst, we utilized 
the spectral data from the "real time" telemetry modes. For these bursts, spectral data 
were available with either 2, 16, or 64 second accumulations. For the longer accumulations, 
spectral studies were not generally feasible. See Murakami et al. (1989) for more information 
concerning the detectors. 
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Before a model can be fitted to the data the background spectrum must be estimated 
and subtracted from the counts. For a majority of the analyzed bursts a linear fit to 
the MRO data in each energy channel provided a reasonable fit to the pre- and post- 
burst data. In a few cases, a quadratic fit was used. Events with large variations in 
the background were rejected from this analysis. For several events, the background was 
estimated from real-time data which was generally available for a longer time period both 
before and after the event. The background subtracted spectra were then fitted using a 
standard minimization technique. In this method the input model spectrum is first 
folded through the detector response functions, thereby providing count predictions for 
each energy channel, these predictions are then compared to the observed counts using a 
statistic. The model parameters are adjusted iteratively until a minimum in x^ is obtained. 
As mentioned previously, we have adopted the spectral model employed by B93 because of 
its relative simplicity and ability to accurately characterize a wide range of spectral continua, 
in addition this choice facilitates direct comparison of our results with those obtained from 
the BATSE bursts. This model has the form 



where A is an overall scale factor, a is the low energy spectral index, /3 is the high energy 
spectral index, and Eq is the exponential cutoff or bend energy. In addition, to facilitate 
comparison with previous calculations we have restricted the ranges of /? to be greater 
than -5.0. Recently, this spectral form has been characterized with the peak of the uFj^ 
distribution (i.e., Ep). The conversion is that Ep = {a + 2)Eq. 

The Ginga GBD was in operation from March 1987 to October 1991. During this 
time ~ 120 7-ray bursts were identified (c/. Ogasaka et al. 1991; & Fenimore et al. 1993). 
We selected from this group a sample of 22 events for which good spectral data were 
available and for which we could be reasonably certain that the burst occurred within the 
forward, tt steradian field of view of the detectors (front-side events). Front-side events are 
easily identified because they have consistent fluxes in the energy range observed by both 
instruments (15 to 25 keV). This sample includes 18 of the events searched for spectral 
lines (Fenimore et al. 1993). Excluded bursts usually show strong rolloffs below 25 keV 
and spectral fits with an absorption component due to aluminum (a principal spacecraft 
material) are consistent with this interpretation. Sky positions for four of the events in 
our sample are known because of simultaneous detections with either BATSE or WATCH. 
Incidence angles into GBD for bursts 910429, 910717, and 910814 (9.3, 34.3, 37.0 degrees, 
respectively) were derived from knowledge of the GBD orientation at the time of the burst 
and the published BATSE positions (Fishman et al. 1994). Burst 900126 was seen by the 
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WATCH experiment and the incidence angle (50 degrees) was derived using the pubhshed 
position (Lund 1992). 

For the remaining events, the incidence angle 9 of the photons into the detectors is 
uncertain (0° equals normal incidence). Since the detector response is a function of this 
angle, the inferred source spectrum and thus peak intensity of these events is also somewhat 
uncertain. We selected 37 degrees for the incidence angle when the angle was unknown. 
This is a typical angle considering that the mechanical support for the window on the PC 
acts as a collimator limiting the field of view to an opening angle of ~ ±60 degrees. We 
have used Monte Carlo simulations to evaluate the impact of the unknown incidence angle. 
For each simulation, we selected a random incidence angle between and 60 degrees. We 
calculated a response matrix and used it with the burst's best fit parameters in Equation (1) 
to generate simulated spectra. Background was added and Poisson statistics applied. We 
then analyzed the simulated data the same way we analyzed the burst data: an estimated 
background was subtracted and the best fit parameters for Equation (1) were found based 
on a response matrix corresponding to 37 degrees. 

The uncertain incidence angle usually only affects a and Eq in Equation (1) so we 
will emphasize the affect of the uncertain angle on those parameters. Figure 1 shows 
typical examples covering interesting values of a (i.e., positive a's: 0.74, 0.22, and 1.67 for 
GB900901, GB870303, and GB900322, respectively). Each cross represents the best fit a 
and Eq from a simulation. We show 100 simulations per burst for clarity although we did 
10"^ simulations. The contour includes 68% of the cases (based on the 10"^ simulations). 
As was found by B93, the observations tend to agree with a range of solutions satisfying 
a oc logE'o. Figure 1 shows the range of uncertainty due to both the unknown incidence 
angle and counting statistics. 

3. X-ray Spectral Characteristics of GINGA Bursts 

For each burst in the sample we selected a time interval AT over which to sum the 
counts in each spectral channel. This selection was based on the need to increase the 
signal to noise in the computed spectrum and the desire to constrain the spectrum near 
the peak of the time profile. Sometimes, AT is limited by the available time resolution 
in the real time data. If the burst has a smaller duration than our resolution, we will 
underestimate its intensity (e.g., GB910814). Spectra can evolve substantially during the 
course of an event (c/. Mazcts et al. 1981; Golcnetskii et al. 1983; Norris et al. 1986; Band 
et al. 1992; Kargatis et al. 1994; Ford et al. 1995; and §4). We attempted to compute a 
spectrum consistent with the most intense portion of each burst's time profile. For real-time 
events we endeavored to compute the spectra from the time ranges which encompassed the 
majority of the burst. We employed all spectral channels in most of the model fits except for 
the lowest PC channel which has extremely small effective area, for a total of 47 channels. 
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To illustrate these spectral fits we show in Figure 2 the best-fit spectrum and indicate 
the contribution of for each channel for each burst analyzed. The burst nomenclature 
specifies the year /month/day on which each event was detected. In the event that more 
than one burst occurred on a given day, a letter is assigned to the burst according to the 
order of occurrence. 

We choose the format of Figure 2 to simultaneously show the best-fit spectra as well 
as how well each spectral channel fits the observations. Note that the ordinate is spectral 
number flux, not counts. Let Oj be the observed net count rate for the i-th channel and 
Mj be the predicted net rate from the model, N{E). Let af be the variance on Oj. The 
solid line in these figures is N{E) and each channel is represented by a point plotted at 
OiN{E)/Mi with a vertical line length of ±aiN{E)/Mi. Thus, for each point, the distance 
from the point to the line divided by the length of the error bar is (Mj — Oj)/(Tj = Xj, that 
is, it gives the contribution of the point to the total x^. The alternative ways to present 
these data would be to plot the residuals Mj — Oj it or to plot Mj and Oj ± Uj. Each 
method emphasizes a different aspect of the fitting: the model photon spectrum by the first, 
the residuals by the second, and the observed counts by the third. It is equally possible 
to visualize the contribution of each channel to i^i three presentations. All three 
methods are only as valid as the assumed N{E). The third method is often characterized 
as being independent of the model but, of course, the only useful information is whether the 
observed counts differ significantly from the predicted counts, and those predicted counts 
depend on the model. One most be careful in all three presentations to avoid interpreting 
deviations of Oj from Mj as implying that the spectrum might be locally different from the 
assumed N{E) because is all three cases, Mj depends on a wide range (i.e., from E to oo) 
of N{E) so is very model dependent. The temptation to view local differences between Oj 
and Mj as evidence that N{E) should be different is probably strongest for our format since 
N{E) is present in the figure. However, all three formats are equivalent on this score: the 
quality of a model fit should be judged solely on the basis of the total value and not on 
how any one point (or group of points) matches (or mismatches) the best-fit model. Our 
format is misleading if one attempts to interpret the deviations as implying that the best-fit 
model should locally change to "agree" with the observations. That is why this format is 
inappropriate for presentations of absorption line candidates or to argue that a model needs 
to be changed. However, in this work, we are showing acceptable fits and we will not argue 
for changes in the best fit model. This format shows the best fit spectral function and how 
well each point agrees with it. As such, it is a good method for visualizing the low energy 
behavior of GRBs. 

The results of the spectral fits for each burst are summarized in Table 1. Here, AT 
specifies the time interval for calculation of the spectrum, Rx/g is the ratio of the emitted 
energy in x-rays relative to the gamma-rays (see §5), Eq, a, and P are the best-fit model 
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parameters, and Xr reduced x for the fit. Each spectral fit was performed with 

43 degrees of freedom, and the quoted uncertainties for (3 are la estimates assuming one 
parameter of interest. The strong coupling in the B93 spectral shape between a and Eq 
means that one should not quote separate confidence regions for a and Eq (see below). 

Previous spectral analyses have been reported for some of the bursts listed here, mostly 
with regard to the presence of lines in their spectra, but also because of interesting X-ray 
properties (c/. Murakami et al. 1988; Fenimore et al. 1988; Wang et al. 1989; Yoshida et al. 
1992; Yoshida et a/. 1989; &; Murakami et al. 1991), here we have restricted our interest to 
the spectral continua of the bursts. In Table 1 we also give references to previous analyses 
of these bursts (last column of Table 1). 

The spectral fits are generally acceptable, with Xr of order unity for most of the bursts. 
In agreement with B93 we find that a range in these model parameter values is required 
to adequately describe 7-ray burst spectra. In Figure 3 we show a graphical comparison 
between the BATSE results taken from B93 and our results for the Ginga sample. Figure 
3a contains the best-fit spectra for the 22 bursts in the Ginga sample, while Figure 3b 
displays the 54 spectra computed from the BATSE data of B93. In each case the spectra 
are plotted only over the nominal band-pass of each instrument (2-400 keV and 20-3000 
keV for Ginga and BATSE respectively). For BATSE, bursts were measured over slightly 
different bandpasses due to different detectors having different gains. See B93 for the actual 
measured range of each burst. The spectra have been normalized to 1.0 photon/keV at 100 
kcV. Between 20 and 400 keV the two samples are in substantial agreement, and span a 
consistent range of spectral shapes. 

Of particular interest is the behavior of the Ginga sample at X-ray energies. About 
40% of the bursts in the sample show a positive spectral number index below 20 keV (i.e., 
a > 0), with the suggestion of rolloff toward lower energies in a few of the bursts {a 
as large as ~ -1-1.5). Unfortunately, the lack of data below 1 keV, and the often weak 
signal below 5-10 keV precludes us from establishing the physical process (photoelectric 
absorption, self-absorption) that may be involved in specific bursts. The remainder of the 
burst spectra continue to increase below 10 keV. Observations of the low-energy asymptote 
can place serious constraints on several GRB models, most notably the synchrotron shock 
model which predicts that a should be between -3/2 and -1/2 (Katz 1994). Crider et al. 
(1997) uses examples from BATSE to argue that some GRB violate these limits during 
some time-resolved samples. We find violations of these limit in the time-integrated events. 

Figure 4 summarizes the low energy behavior of GRBs as a scatter plot of the bursts 
a's and Eq^s. In Figure 4a, we show the 68% confidence regions for 20 of the bursts in 
Table 1. These confidence regions were calculated the same way as in Figure 1. For two 
bursts we only show the best fit parameters as solid squares because a simple power law 
could nearly fit the entire spectra. Note for GB881009 that a is nearly equal to [3 and for 
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GB880205, the bend energy is well above our energy range. Although we give the formal 

best fit parameters in Table 1, effectively, Eq was undetermined or not necessary for the 
fit in those two bursts. We repeated the analysis of Figure 4a in Figure 4b except we did 
each simulation at the angle for which we analyzed the burst (37 degrees, in most cases). 
Thus, Figure 4b indicates the size of the confidence regions if we knew the incidence angle. 
By comparing Figure 4a and Figure 4b, we see that the lack of knowledge of the incidence 
angle into Ginga does not introduce much more uncertainty than the counting statistics. 
On average, because of the uncertain incidence angle, the confidence region for Eq is 22% 
larger and the confidence region for a is larger by 0.06. 

In Figure 4c we combine the results from this paper and B93. The open squares are the 
BATSE results from B93 and the solid squares are the Ginga results presented here. Many 
of the Ginga points lie within the range found by BATSE. However, the lowest Eq found 
by BATSE was 14 keV (set, of course, by the lowest energy observed by BATSE). Ginga 
extends Eq values down to 2 keV. BATSE had a small fraction (15%) of events with a > 
whereas Ginga has 40% of events with a > 0. In general this is because there is a correlation 
between a and Eq such that the lower energy range of Ginga samples a parameter space 
with more events with a > 0. For the 76 points, the Pearson's r coefficient is -0.62. The 
formal significance is about 4cr although that ignores the complicated error bars that are 
caused by the fact that the observations tend to agree with a range of a — Eq. However, the 
existence of a correlation seems reasonable: there are virtually no events seen by BATSE 
at large a, large Eq and few low a, low Eq events seen by Ginga. The error bars (Fig. 4a) 
are large on the Ginga points, but not large enough to indicate that they all should fall 
where the average BATSE events occurred. We note that GB900126 is an event for which 
we know the incidence angle. This event is important to our conclusions because it has one 
of the largest a's and there is no uncertainty due to an unknown incidence angle. In Figure 
4c we show a representative confidence regions for one BATSE event from B93, burst 451. 
We note that the confidence regions quoted in Table 4 of B93 were found separately for a 
and £^0) ignoring the coupling, so are much smaller than the real confidence regions and 
cannot be compared to ours. 

In our sample of bursts, we find that a can be both less than zero and greater than 
zero. Negative a's arc often seen in time-integrated BATSE spectra. Positive a's where the 
spectrum rolls over at low energies are usually only seen in time resolved BATSE spectra 
(Crider et al. 1997). 

The Ginga trigger range (50 to 400 keV) was virtually the same as BATSE's. Thus, we 
do not think we are sampling a different population of bursts, yet we get a different range 
of fit parameters. The lack of events with Eq^s between 6 and 20 keV cannot be used to 
support two populations because we do not have enough events. One possible explanation 
might be that GRBs have two break energies, one often in the 50 to 500 keV range and the 
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other near 5 keV. Both BATSE and Ginga fit with only a single break energy so BATSE 
tends to find breaks near the center of its energy range and we tend to find breaks in our 
energy range. Without good high energy observations of bursts with low Eq, it is difficult 
to know whether they also have a high energy bend. 

Preece et al (1996) utilized a low energy discriminator channel and detected emission 
in excess of what would be expected from a fit at higher energy. Preece et al (1996) report 
excesses in 15% of the investigated BATSE bursts. One of our bursts, GB880205, shows a 
clear strong excess at low energy and two other Ginga bursts, GB880830 and GB910418, 
probably also shows an excess. For GR910418, Preece et al (1996) also reported an excess. 
Prom the Preece et al (1996) result, wc would expect about 3 of our bursts to show an 
excess so we are consistent with the BATSE result. 

4. Spectral Evolution of GINGA Bursts 

It was recognized shortly after the discovery of the 7-ray bursts that their spectra 
show significant variation within a single event (Wheaton et al. 1973; Mctzgcr et al. 1974; 
Teegarden & Cline 1980; &: Yoshida et al. 1989). Analysis of Solar Maximum Mission, 
and BATSE data indicate that burst spectra tend to evolve from hard to soft, both for 
individual peaks in the temporal profile, as well as over the entire course of the burst, 
though counter-examples to this general trend can be found (c/. Norris et al. 1986; and 
Band et al. 1992, Band, 1997, Crider et al. 1997). Previous analyses of the Ginga data 
identified several bursts with soft X-ray preactivity as well as X-ray tails which generally 
showed softer spectra than the most intense portion of the event (Murakami et al. 1991; 
and Murakami et al. 1992b). 

Several of the Ginga bursts which wc have analyzed were of sufficient duration and 
intensity to permit spectral evolution studies. Here we summarize our results for two 
events, bursts 901001 and 890929. Incidence angles into the detectors for these events are 
not known, however, when comparing spectra computed from different time intervals within 
the same burst, the choice of incidence angle is not important, since this angle did not vary 
during the event. We therefore used an angle of 37° for these calculations as well. For each 
burst we subdivided the time profiles into several temporal regions and then computed the 
best-fit spectrum for each region. We used the same GRB model (eq. [1]) as in the previous 
calculations. 

In Figure 5a we show the time profiles for burst 901001 in five different energy ranges. 
The dotted vertical lines delineate the six temporal regions in which spectra were computed. 
Also shown is a linear estimate of the background computed from the pre- and post-burst 
signal. Note that the low energy. X-ray emission persists longer than the high energy 
emission, suggestive of overall softening of the burst spectrum {of. Yoshida et al. 1992). 
In Figure 5b we show the best-fit spectra computed for each indicated time period. We 
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have again normalized each spectrum to 1.0 at 100 keV. The Figure legend identifies each 
spectrum with the time interval for which it was computed. For completeness, the best-fit 
model parameters computed from each time interval are also given in Table 2. Comparison 
of the spectra computed from the first two time intervals suggests possible initial hardening 
during the rise to peak intensity. After the peak there is strong evidence for softening 
of the spectrum. In fact, for intervals and Tq there is still a significant signal in the 
PC, however, there is no evidence of corresponding emission above 100 keV in the SC. To 
further investigate the relationship between intensity and hardness we have computed Ep 
and flux F (between 2 and 400 keV) for each time interval. The resulting values are plotted 
in the inset panel of Figure 5b. The la uncertainties were estimated from Monte Carlo 
simulations using the best-fit spectral parameters derived from each time interval. The 
direction of temporal evolution is along the solid line, from upper right to lower left. Thus, 
there is initial hardening during the rise to peak of the burst, followed by softening for the 
remainder of the event. In Figures 6a, 6b and Table 2 we show the corresponding results 
for burst 890929. There is also evidence for an X-ray tail in this event, and the softening 
of the spectrum during the course of this burst is also quite apparent, though no initial 
hardening during the rise to peak is evident. For these events the X-ray signal also softens 
in addition to the higher energy 7-rays. 

The hardness-intensity plots for these two bursts suggest a correlation between hardness 
and intensity similar to that reported by Golcnetskii et al. (1983) , Kargatis et al. (1994), 
and Bhat et al. (1994). Linear fits to the logEp versus logF data (i.e. logEp = alogF + b) 
for bursts 901001 and 890929 give a « 0.3. The results of our linear fits are shown in Table 
3. 

5. X-ray Emission Relative to Gamma-ray Emission 

Early measurements of the x-rays associated with gamma-rays were fortuitous obser- 
vations by collimated x-ray detectors that just happen to catch a GRB in their field of 
view. From a few events it appeared that the amount of energy in x-rays was only a few 
percent (Laros et al. 1984) confirming that GRBs were, indeed, a gamma-ray phenomena. 
The Ginga experiment was designed with a wide field of view to detect a sufficient number 
of events to determine the range of x-ray characteristics. Early reports from Ginga events 
indicated that sometimes a much larger fraction of the emitted energy was contained in the 
x-rays. For example, by comparing the signal in the proportional counter (roughly 2 to 25 
kev) to that of the scintillator (roughly 15 to 400 kev), we reported an x-ray to gamma-ray 
emission ratio up to ~46% (Yoshida et al. 1989). Such a ratio depends on the bandpass for 
which it is evaluated. One can use the parameters from Table 1 in Equation (1) to determine 
the ratio of emission for arbitrary bandpasses. This would be useful, for example, to predict 
the range of emission that might be seen in instruments such as ROSAT. Here, we find the 
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ratio of emission for a typical x-ray bandpass (2 to 10 keV) compared to the BATSE energy 
range (50 to 300 keV). The ratio is defined to be R^jg = EN{E)dE / ^^^^ EN{E)dE. 
This ratio is listed in the third column of Table 1. Figure 7 presents the distribution for the 
events analyzed in this paper. Although the ratio is often a few percent, for some events 
the ratio is near (or larger than) unity. Some GRBs actually have more energy in the x- 
ray bandpass than the gamma-ray bandpass. The simply average of the 22 values is 24%. 
This large value arises because of the few events with nearly equal energy in the x-ray and 
gamma-ray bandpass. However, even the logarithmic average is 7%. 

6. Discussion and Summary 

Several of the results reported here are of particular interest. There is some concern in 
cosmological models of GRBs that the bend energy occurs only over a small range of values 
implying that the bulk Lorentz factor only varies a little from burst to burst (Brainerd, 
1994, 1997). The Ginga bursts extends the range of observed bend energies from above 500 
keV (with BATSE) to 3 keV (Ginga). This Ginga result implies that, perhaps, we have 
not yet observed the lower limit on bend energies and, thus, the range of dynamic range of 
bulk Lorentz factors has not been limited. Regardless, any successful model for the burst 
sources must be able to explain the range of bend energies. 

The spectral evolution observed in bursts 901001 and 890929 support the consensus 
that hard to soft spectral evolution is a prevalent feature of the burst mechanism. Both 
Ginga bursts analysed here would seem to fit this scenario. Initial soft to hard evolution 
has been previously reported for burst 900126 (Murakami et al. 1991), and is also suggested 
by our analysis of burst 901001. Perhaps the timescale for the energization process is rapid 
enough in some bursts to preclude observation of the initial hardening of the spectrum. 

Several future space missions have been planned with the goal of detecting bursts at 
soft X-ray energies. For example, the X-ray spectral distribution of 7-ray bursts has conse- 
quences for the High Energy Transient Experiment (HETE) mission (c/. Ricker et al. 1992) 
which has as one of its main goals the identification of burst counterparts derived from 
localization of bursts using X-ray observations. The initial localization of sources by HETE 
depends on the detection of bursts with a coded aperture X-ray instrument operating from 
2-25 keV. Estimates of the number of bursts expected to be seen per year can be made 
based on the range of x-ray spectra presented here. 

We have investigated the X-ray spectral characteristics of a sample of 22 bright 7-ray 
bursts detected by the GBD which flew aboard Ginga. In the energy band where they can be 
compared, the range of spectral shapes which describe the Ginga sample are consistent with 
those derived from analyses of bright BATSE bursts. Consistent with Preece et al. 1996, 
we find at least one event (GB880205) that has an x-ray excess relative to the GRB spectra 
from Equation 1. Moreover, the Ginga spectra extend down to 2 keV, and thus provide 
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a useful characterization of burst spectra over an energy range at which burst spectra have 
not often been measured. In addition, we have investigated spectral evolution in two longer 
duration events. We found evidence for spectral softening during both of these events, along 
with a correlation between hardness and intensity. 

This work was carried out under the auspices of the US Department of Energy. It is a 
pleasure to acknowledge helpful discussions with Jean in 't Zand and Richard Epstein, and 
Tony Crider. We also thank David Band for extensive discussions on the error bars for a 
and ^ well as many comments on the manuscript. Tony Crider kindly helped with the 
production of Figures 5 and 6. 
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Figure Captions 

Figure 1-Simulations for three bursts to determine the effects of the unknown incidence 
angle and the counting statistics. Each cross is a simulation for a burst located at an 
angle randomly selected between and 60 degrees but analyzed with a response matrix 
corresponding to 37 degrees, the response matrix used to analyze the Ginga spectra. One 
hundred simulations are shown for clarity but 10^ were used. The contour includes 68% of 
the events. 

Figure 2-Best fit spectrum for each burst in our sample. The vertical error bars are ilcr, 
and the distance from the photon spectrum to the center of the horizontal bar represents 
the residual in units of la. Thus, the position of the points represent the contribution of 
that data point to the of the fit (see text) . The horizontal bars represent the energy loss 
bins. 

Figure 3a-The best-fit spectra for the 22 Ginga bursts in our sample. Each curve has been 
normalized to 1.0 photons/keV at 100 keV. The spectra are plotted from 2.0 to 400.0 keV, 
the approximate range of sensitivity of the Ginga GBD. 

Figure 3b-The best-fit spectra for 54 BATSE bursts. The best-fit parameters are taken 
from B93. We have plotted the spectra from 20.0 keV to 3 Mev, the approximate band pass 
for the BATSE spectral fits. 

Figure 4a-Simulations to determine the 68% confidence regions of a and Eq. The confidence 
regions were found in the same manner as in Figure 1 and include both the effects of the 
unknown incidence angle and the counting statistics. The nature of the model (eq. [1]) is 
such that a range of parameters satisfying a cx log Eq can be consistent with the data. For 
two events, Eq was unnecessary or undetermined and they are shown as solid squares. 

Figure 4b-Confidencc regions reflecting the uncertainty due just to the counting statistics. 
For these simulations it was assumed that the incidence angle was known and the same as 
the analysis angle (37 degrees, for most bursts). By comparing Figure 4a and Figure 4b, 
we see that the effects of the uncertain incidence angle is minor compared to the effects of 
the counting statistics. The uncertain angle increases, on average, log(^o) by 0.09 and a 
by 0.06. 

Figure 4c-Low energy slopes and bend energies for Ginga and BATSE events. The open 
squares are 54 BATSE events from B93 and the solid squares are the 22 Ginga events 
reported in this paper. The Ginga events extend the observed range of break energies to 
lower values and reveal a correlation (Pearson's r = —0.62) between the low energy slopes 
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(a) and the bend energy {Eq). For comparisons with Figure 4a, we show the confidence 
region for one BATSE event (number 451). 

Figure 5a-Time history of burst 901001 recorded in five energy bands, two from the PC 
and three from the SC. The dotted vertical fines denote the time intervals in which spectra 
were computed. The dashed line is a linear estimate of the background derived from a fit to 
pre- and post-burst data. To displace each profile vertically, a constant value of 600 counts 
has been added successively to each profile. 

Figure 5b-Spectra computed from the six time intervals of burst 901001 delineated in Figure 
5a. The legend indicates the time interval for which each spectrum was computed. The 
inset Figure shows the hardness-intensity evolution of this event. The direction of temporal 
evolution is along the solid line, from upper right to lower left. 

Figure 6a-Time history of burst 890929 recorded in five energy bands, two from the PC 
and three from the SC. The dotted vertical line denote the time intervals in which spectra 
were computed. The dashed line is a linear estimate of the background derived from a fit to 
pre- and post-burst data. To displace each profile vertically, a constant value of 600 counts 
has been added successively to each profile. 

Figure 6b-Spectra computed from the six time intervals of burst 890929 delineated in Figure 
6a. The legend indicates the time interval for which each spectrum was computed. The 
inset Figure shows the hardness-intensity evolution of this event. The direction of temporal 
evolution is along the solid line, from upper right to lower left. 

Figure 7-The distribution of the ratio of the energy emitted in x-rays relative to that emitted 
in gamma-rays. The x-ray bandpass is defined to be from 2 to 10 kev and the gamma-ray 
bandpass is the BATSE range of 50 to 300 keV. Note there are some examples of equal 
energy in the x-rays and the gamma-rays. The average ratio of the energy in the x-rays to 
the energy in the gamma-rays is 24%. 
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^ Year /month/day of event 

^ Accumulation time of spectrum in sec 

ratio of the emitted energy in 2 to 10 kev relative to 50 to 300 kev (see § 5) 
^ Break energy for Equation 1 (in kcV), uncertainties are shown in Figure 4a 
® Low energy spectral index, uncertainties are shown in Figure 4a 

High energy spectral index 
S Reduced for best fit spectrum (43 dof) 

^ References to previous work 1: Murakami et al. 1988, 2: Fenimore et al. 1988, 3: Yoshida 
et al. 1989, 4: Wang et al. 1989, 5: Murakami et al. 1991, 6: Graziani et al. 1992, 7: Lund 
1992, 8: Yoshida et al. 1992, 9: Murakami et al. 1992a, 10: Murakami et al. 1992b, 11: 
Graziani et al. 1993, 12: Yoshida & Murakami: 1994, 13: Fishman et al. 1994, 14: Schaefer 
et al. 1994 15: Fenimore et al. 1993, 16: Band et al. 1993, Preece et al. 1996. 
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Table 2 

Spectral Evolution Parameters* 
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* See § 2 for a discussion of the parameters. 



Table 3 

Hardness-Intensity Correlations* 

Bursts a aa b ai, P 

901001 0.276 0.035 3.79 0.208 0.05 

890929 0.471 0.042 4.90 0.228 0.13 



* Fit of the form log Ep = a log F + b 
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